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♦p^l We have identified 335 galaxy cluster and group candidates, 106 of which are 

O I at 2 > 1, using a 4.5/im selected sample of objects from a 7.25 deg^ region in the 

^ [ Spitzer Infrared Array Camera (IRAC) Shallow Survey. Clusters were identified 

' as 3-dimensional overdensities using a wavelet algorithm, based on photometric 
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redshift probability distributions derived from IRAC and NOAO Deep Wide- 
Field Survey data. We estimate only ~ 10% of the detections are spurious. To 
date 12 of the z > 1 candidates have been confirmed spectroscopically, at redshifts 
from 1.06 to 1.41. Velocity dispersions of ~ 750 km s~^ for two of these argue 
for total cluster masses well above 10^'^ Mq, as does the mass estimated from the 
rest frame near infrared stellar luminosity. Although not selected to contain a 
red sequence, some evidence for red sequences is present in the spectroscopically 
confirmed clusters, and brighter galaxies are systematically redder than the mean 
galaxy color in clusters at all redshifts. The mean I — [3.6] color for cluster 
galaxies up to 2; ~ 1 is well matched by a passively evolving model in which stars 
are formed in a 0.1 Gyr burst starting at redshift Zf = 3. At 2; > 1, a wider range 
of formation histories is needed, but higher formation redshifts (i.e. Zf > 3) are 
favored for most clusters. 

Subject headings: galaxies: clusters: general — infrared: surveys 



Introduction 



As the most massive gravitationally bound systems in the Universe, the rate of emer- 
gence of galaxy clusters since the Big Bang might be expected to be among the most straight- 
forward predictions of cosmological models. Yet despite the advent of the era of precision 
cosmology ushered in by observations of SNe and the cosmic microwave background (CMB), 
significant uncertainty remains in the expected numbers of galaxy clusters at z > 1. The 
CMB temperature anisotropies on scales corresponding to clusters are not accurately known, 
leading to a range of values for erg, the rms matter fluctuation in a sphere of radius 8h~ ^ 
Mpc at ;z = 0. Est imates of a^ vary s ignificantly, i ncluding e.g. 0.67"*"^ j^ (JGladders et al. 



20071). 0.76 ± 0.05 flSpergel et all 120071). 0.80 ± 0.1 fJHetterscheidt et al 



I 



200711. 0.85 ± 0.06 



(iHoekstra et al. 120061). 0.90 ± 1 f|Spergel et al.l l2003h. 0.92 ± 0.03 flHoekstra et all 120021 ) . 
and 0.98 ± 0.1 (JBahcall Sz Bodell2003l ). The range erg = 0.7 to 1 corresponds to a variation 



of a factor of nearly 20 in the predicted numbers of 2; > 1 clusters with Mtot > 10 M© (e.g 



Sheth fc Tormerull999l ). Removing this uncertainty is a major goal of upcoming Su nyaev- 
Zeldovich cluster surveys such a s the SZA (Sunyaey -Zeldovich Array: iLoh et al.ll2005l ). AMI 
(Arcminute Mic rokelvin Imager: iKneissl et al.ll2001l). ACT fAtacain a Cosmology Telescope; 
Kosowskvl[2003l ). and SPT (South Pole Telescope: IRuM et al.l[20oi ). 

Since by definitioiu galaxy clusters contain an unusually high density of galaxies, they 



^There is no generally agreed upon definition for galaxy clusters: see e.g. lAbelll (|l958l ). iPostman et al 
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provide an efficient means of observing substantial numbers of galaxies at a common distance, 
offering the hope of constructing the analog of the Hertzsprung- Russell (i.e. color-magnitude) 
diagram for galaxy evolution. Indeed, studies of the relationship between color and mag- 
nitude indicate that clu sters are the habitat of galaxies with the oldest and most ma ssive 
stellar populations (e.g. [Stanford. Eisenhardt. fc Dickinsonlll998l : iBlakeslee et al.ll2003l ). ob- 
jects reasonably free of the complications associated with starbursts and dust. These studies 
are consistent with an extremely simple formation history for cluster galaxies, in which their 
stars are formed in a short burst at high redshift, and they evolve quiescently thereafter (we 
use the term "red spike model" in referring to this scenario - see Figure [J). Studies of the 



near-IR luminosity fuiictions of cluster galax i es reinforce this pictu re (e.g. Ide Propris et al. 
I999I : be Propris et al.lboOTi : IXoft et al.lEooi : IstrazzuUo eraPbood ). With their large look- 
back times, therefore, high redshift galaxy clusters also provide an observational pillar for 
our understanding of the formation and evolution of galaxies. 

Obtaining substantial samples of galaxy clusters at 2; > 1 has proved challenging, largely 
because such objects are difficult to detect using only optical data. Due to their greatly en- 
hanced rate of star formation by 2; ~ 1, the UV emission from modest sized field galaxies 
overwhelms that from the intrinsically red spectra of quiesecen t, early type galaxies prefer- 
entially found in clusters. The Red Sequence Cluster Survey ( IGladders fc Yed l2000l . l2005l ) 
uses the observed color-magnitude relationship in cluster galaxies to improve the contrast 
and has proven highly efficient to z ~ 1, but the optical colors of the red sequence become 
increasingly degene rate at higher redshifts, as they no longer span the rest 4000A break. 
Wilson et al.l (120061 ) describe a program to extend the red sequence technique to higher red- 
shift using Spitzer data, but it is also important to test for the existence of red sequences in 
z > 1 clusters rather than preselecting for them, if possible. 

The contrast of high redshift clusters over the field improves at longer wavelengths 
(Figured]), but the contrast against atmospheric emission declines, and until recently the 
relatively small formats of infrared detect o r arra ys made surveying sufficient z > 1 volume 
a formidable undertaking. [Stanford et al.l (119971 ) reported the discovery of a cluster at z = 
1.27 in a 100 square arcminute survey to KsiVega) = 20 (lOa). But this survey required 
approximately 2 hours of exposure in both J and Kg per position, and 30 allocated nights 
of KPNO 4m time to complete. With estimates for the surface density of IO^^Mq clusters 
at 2; > 1 in the range 0.2 — 4 per square degree (ISheth fc TormenI 119991 ). the discovery was 
in hindsight fortuitous. 



(|l996l ). and lRosati et al.l (|2002l ). We define our criteria for candidate galaxy clusters and groups in ^ and 
for spectroscopically confirmed candidates in Sj6l 
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Fig. 1. — L* for cluster galaxies vs. z, in the observed By/ (long dashed ), / (short dashed), 
Kf, (d otted), and IRAC 4.5/im (solid) bands. The curves are based on a lBruzual fc Chariot 
( 120031 ) model where stars are formed in a 0.1 Gyr burst beginning at 2/ = 3 in a A = 0.7, 
O m = 0.3, h = 0.7 c osmology, which fits the observed L* in galaxy clusters to 2; > 1 



( Ide Propris et al.l Il999h . This model is referred to as the "red spike" model in the text. 
Horizontal lines show the 5a limits of the IRAC Shallow Survey, FLAMEX survey, and 
NDWFS in 5" diameter apertures. Note a 5" aperture is larger than optimum for detection 
at B\y and /, and the 50% completeness limits in these bands are about 2 mag fainter. 
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Such considerations motivated a different approach, where extended sources in deep X- 
ray surveys lacking prominent optical counterparts were targeted fo r IR foUowup. This 
technique yielded confirmed clusters ai z = 1.10, 1.23, and 1.26 ([Stanford et al.l 12002 



Rosati et al.lll999l . |200J). With the arrival of XMM, X-ray surveys offer renewed pro mise, 
leadir ig recently to the ide ntification of galaxy clusters a.t z = 1.39 (JMuUis et al.ll2005l ) and 
1.45 (jStanford et al.l 120061 ). With exposure times > 20 ksec and a 30 arcmin field of view, 
a discovery rate of approximately 30 hours per candidate z > 1 cluster above lO^^M© is 
expected (assuming one such cluster per square degree, which corresponds to ag = 0.83). 

Searches fo r clusters around rad i o galaxies have yielde d prot oclusters with redshifts 
as high as 4.1 (IPentericci et al.l l2000l : IVenemans et al.l 120021 . |2005| ) and possibly even 5.2 
( lOverzier et al.ll2006l ). The very large redshifts of these systems enable powerful inferences 
to be drawn regarding the formation of cluster galaxies, but they are less useful as probes 
of the cosmological growth of structure. Lyman Break Galaxy surveys with intensive fol- 
lowup spectroscopy o n the Keck telescopes ha ve al so identified highly ove rdense structures 
atz = 2.30 and 3.09 Jsteidel et al.lll998l . l2005h . and lOuchi fc SXDSl J2007h discuss a ^ = 5.7 
structure identified via Ly-a emission in a narrow band imaging survey. 

Recent advancements in IR detector array formats have renewed interest in ground- 
based IR surveys. In one example of the state of the art, Elston, Gonzalez et al. (2006) use 
the 2048 x 2048 pixel FLAMINGOS camera to map 4 deg^ to a 50% completeness limit of 
Kg = 19.2 (Vega). With 2 hour exposures on the KPNO 2.1m each covering 1/lOth deg^, 
this leads to an expected discovery rate for high redshift clusters per useful hour of observing 
which is similar to XMM. The UKIDSS Ultra Deep Survey provides ano ther recent example, 
findin g 13 cluster candidates with 0.6 < z < 1.4 in a 0.5 deg^ s urvev (jvan Breukelen et al. 
20061). one of which h as 4 spectroscopic redshifts at z = 0.93 (JYamada et al.ll2005l ). and 



Zatloukal et al.l (120071 ) find 1 2 candidates with 1.23 < z < 1.55 in a 0.66 deg^ //-band 
survey in the COSMOS field. iMcCarthy et al.l (120071 ) present a system with a high density 
of galaxies with red optical to near-IR colors surrounding a galaxy at 2; = 1.51, identified in 
the 120 arcmin^ Gemini Deep Deep Survey. Candidates drawn from surveys of less than a 
square degree are unlikely to include many rich clusters, however. 



With the launch of the Spitzer Space Telescope in 2003 (I Werner et al.l 120041), sensitive 



infrare d arrays free from foreground thermal emission were put into operation ( 



Fazio et al. 



2004bl ). A major scientific driv er for the Spitzer Infrared Array Camera (IRAC) Shallow 
Survey (JEisenhardt et al.l |200J) was the detection of 2; > 1 galaxy clusters. The IRAC 
Shallow Survey uses 90 second exposures per position and covers 8.5 deg^, leading to an 
expected discovery rate of < 8 hours per z > 1 cluster. Here we present results from the 
IRAC Shallow Survey cluster search, finding 106 cluster and group candidates at 2; > 1, of 
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which we estimate only ~ 10% are spurious. 

A surface density of over 10 systems per square degree at z > 1 is higher than expected 
for bound systems with masses above 10^"^ Mq for the current range of plausible Ug estimates. 
While we present evidence that at least two of the z > 1 clusters have masses well above 
1O^^M0, it is likely that our sample includes systems with masses below IO^^Mq (i.e. groups), 
and perhaps some unbound filaments viewed end-on. In the remainder of this paper, for 
brevity the terms "clusters" and "candidates" are used to refer to all such objects which 
meet our selection criteria, unless otherwise stated. 

This paper describes how the cluster sample was identified, and some of the overall 
photometric properties of the sample. We also provide new spectroscopic evidence support- 
ing nine of these clusters, from z = 1.057 to 1.373. Spectroscopic evidence in support of 
IRAC Shallow Sur v ey selected clus t ers a t z = 1.112, 1.243, a nd 1. 413 was presented in 
Elston et al.l (120061 ) . iBrodwin et al.l (120061 ) . and [Stanford et al.l ( 120051 ) respectively, and we 
provi de additional pr e vious ly unpublished spectroscopy on those clusters here, for complete- 



ness. 



Brodwin et al.l (120071 ) discuss the clustering of the clusters, and Galametz et al. (in 



preparation) report on AGN incidence vs cluster-centric distance. FoUowup imaging with 
HST (GO 10496, Perlmutter; 10836, Stanford; and 11002, Eisenhardt) and Spitzer (GO 
30950, Eisenhardt) is underway, and future papers will examine the scatter in the color- 
magnitude relation as a function of morphological type, the dependence of cluster galaxy 
size on redshift, starburst activity in clusters vs. redshift, and the dependence of mean 
galaxy properties on surface density. 

A cosmology with Hq = 7 km s~^, fim = 0.3 , A = 0.7 is assumed, and magnitudes are 
on the Vega system (defined in lReach et al.ll2005l . for IRAC). At 2; = 1 — 2, this means that 
one arcminute corresponds to a physical scale of 480 — 502 kpc, peaking at 508 kpc at z = 1.6. 
Unless otherwise specified, physical (rather than co- moving) scales are used throughout. 



2. Data 



2.1. IRAC Shallow Survey 



The IRAC Shallow Survey (JEisenhardt et al.ll2004j ) was designed to maximize the num- 
ber of reliable sources detected per unit time and to c over sufficient area to det ect significant 
numbers oi z > 1 galaxy clusters. As explained in lEisenhardt et al.l (120041 ). a 30 second 
exposure time per pointing is close to optimum for maximizing source detections, and reli- 
ability was obtained by requiring three independent exposures separated by hours at each 
position. The survey covers ~ 8.5 deg^ and reaches an aperture-corrected 5a depth of ~ 
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19.1 and 18.3 mag (Vega) at 3.6 and 4.5/ini in 3" diameter apertures. It is a remarkable 
fact that 90 seconds of combined exposure with IRAC on the 85 cm Spitzer Space Telescope 
provides sufficient sensitivity to detect evolving L* galaxies to z = 2 (Figure 1). 



2.2. NOAO Deep Wide-Field Survey 

The s urvey was carried out in the Bootes region of the NOAO Deep Wide-Field Survey 
(NDWFS; iJannuzi &: Deylll999l ) to allow photometric redshifts to be derived using the deep 
optical imaging available for this field. The NDWFS reaches 5a point-source depths in the 
B]v,R, and / bands of ~ 27.1, 26.1, and 25.4 respectively. Typical exposure times per 
position with the Mosaic- 1 camera on the KPNO Mayall 4-m telescope were 1-2 hours in 
Bw, 1-2 hours in R, and 2-4 hours in /, and the seeing ranged from 0.7 to 1.5". The data 
acquisition, reduction, and catalog generation are discussed in detail by B. Jannuzi et al. (in 
preparation) and A. Dey et al. (in preparation). This paper uses the NDWFS third data 
release (DR3) images and SExtractor catalogs which can be obtained through the NOAO 
data archive). 



2.3. AGN and Galaxy Evolution Survey, FLAMINGOS Extragalactic Survey, 

and Other Surveys 

The AGN and Galaxy Evolution Survey (AGES, C. Kochanek et al. in preparation) 
provides spectroscopic redshifts for ^ 17, 000 objects (using the version 2.0 catalog) in the 
IRAC Shallow Survey. AGES is highly complete for sources brighter than 15.7 mag at 4.5/im, 
and also for sources brighter than / = 18.5, with many redshifts for sources up to / = 20, 
enabling excellent assessment of the photometric redshifts to z ~ 0.5 (Figure 1). 

Deep near infrared imaging from the FLAMINGOS Extragalactic Survey (FLAMEX) 
is available for half the Bootes region (JElston et al.ll2006l ). and was used for deriving a prior 
on the redshift likelihood functions (see @. 



Imaging of the Bootes NDWFS field has also been obtained i n the radio (Ide Vries et al. 



2OO2I ). at 24, 70, and 160/im with the MIPS instrument on Spitzer ( jHouck et al. 



2005|),inthe 



z— band (jCooll 120071 ). in the UV with GALEX, and in X-rays to a depth o f 5 ksec with the 
ACIS instrument on the Chandra X-ray Observatory (JMurray et al.ll2005l ). but these data 
are not used in this paper. 



^http://www. noao.edu/noao/noaodeep 



3. Sample 

Object selection for the cluster search was carried out in the 4.5/im band, because the 
negative K-correction as the rest-frame l.Qfim peak shifts into this band leads to a flux which 
is nearly independent of redshift fo r 0.7 < z < 2 (Figure 1). Object detection and photometry 



was carried out using SExtractor (JBertin &: Arnoutslll996l ) in double-image mode, allowing 



matched aperture photometry in the other IRAC bands. While smaller apertures maximize 
the depth of the survey, Monte Carlo simulations showed that 5" diameter aperture mag- 
nitudes a re necessary to prov ide sufficiently reliable color measurements and photometric 



redshifts (JBrodwin et al.ll2006l ). This flux limit (5cr at 4.5/im in 5") corresponds to 13.3 /xJy, 



or a Vega-based magnitude of 17.8. 



3.1. Matching to the NDWFS Catalog 

The NDWFS catalogs were also generated using SExtractor, but run in single-image 
mode in each band. Detections in the different optical bands and between the optical and 
IRAC catalogs were matched if the centroids were within 1" of each other, using the closest 
optical source if more than one satisfied this criterion. For very extended objects (generally 
at z < 0.2), detections in the different bands were matched if the centroids were within 
an ellipse defined using the second order moments of the light distribution of the object 



fJBrown et all 120051 ) 



3.2. BvyRI[3.6] Flux Limits and Photometric Errors 

The NDWFS data were taken over several years in variable conditions, and therefore the 
photometric depths vary somewhat from pointing to pointing. Average 50% completeness 
limits for the Bw, R, and /-bands were 26.7, 25.6, and 25.0 respectively. Average 5a fiux 
(magnitude) limits in a 5"aperture (which is significantly larger than the optimum detection 
aperture for these data) were measured via Monte Carlo simulations to be 0.45 /iJy (24.8 



mag) in Bw, 1-03 /iJy (23.7 mag) in R, and 1.45 /iJy (23.1 mag) in the /-band (JBrodwin et al. 



20061 ) . In the space-based IRAC Shallow survey the depth is more uniform, and in the 3.6/im 



band it is 10.0 /xJy (18.6 mag), also derived from a Monte Carlo simulation. 

Although the majority of objects in the 4.5/im-selected catalog are well detected at 
shorter wavelengths, this is not generally the case for z > 1 red ellipticals. Objects at these 
redshifts are often quite faint in the optical as the 4000A break is longward of the /-band 
(see Figure 1). Where sources were observed but not detected, the fiux was taken to be zero 
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and a Monte Carlo la error was adopted. This approach is optimal for photometric redshift 
fitting, where the non-detection provides important contraints on the galaxy spectral energy 
distribution (SED). 



3.3. Sample for Photometric Redshift Estimation and Cluster Search 

The photometry used for photometric redshift estimation (§!]) consists of By[rRI[3.6] [4.5] 
data with the Monte Carlo photometric errors and limits noted above. The 5.8/im and 8.0/im 
bands were not used because they do not have the sensitivity to detect z > 1 cluster L* 
galaxies in the IRAC Shallow Survey. While the IRAC Shallow Survey covers 8.5 deg^ in 
each band, the overlap area observed in both the 3.6 and 4.5/im IRAC bands is 8.0 deg^. 
All of this falls within the optical NDWFS Bootes region. However, due to haloes around 
bright objects, shorter observation times for regions at the edges of individual mosaic camera 
pointings, and residual CCD defects, 0.75 deg^ may have lower quality optical photometry 
and hence photometric redshifts. Hence these regions are also excluded from the sample. 

Finally, 14,044 stars were removed from the catalog using the SExtractor stellarity 
in dex in the bes t seeing optical data for each location. Corn parison to the star count model 
of lArendt et al.l (Il998l ). as tabulated in iFazio et al.l (l2004a[ ). suggests that ~ 80% of stars 
were identified via this approach, leaving ~ 2% of the 4.5yum sample as unrecognized stars. 



In summary, photometric redshifts were estimated for a total of 175,431 objects brighter 
than 13.3 /iJy {5a) at 4.5/xm in 5" in a 7.25 deg^ region. 



4. Photometric Redshifts 



A full description of the photometric redshift methodology is given in iBrodwin et al. 

(120061 ). A summary of those aspects most relevant to cluster detection is provided here. 



Photometric redshifts were computed using an empirical template-fitting algorithm which 
linearly interpolates between the four Cole r nan, W u, & Weedman (1980) SEDs (E, Sbc, Scd, 
and Im), augmented by the iKinney et al.l (119961 ) SB3 ar id SB2 starburst te r nplate s. These 
SEDs were extended to the far-UV and near-IR using iBruzual fc CharlotI (120031 ) models. 
These stellar photospheric models do not include emission from dust, in particular the PAH 
features which dominate the 3 < A < 12yum portion of the spectrum in starforming galaxies, 
but at z > 1 the photometry used does not sample these rest wavelengths. 

In addition to the large AGES (C. Kochanek et al. in preparation) survey, there are 
~ 500 spectroscopic redshifts extending to z ~ 1.5 gleaned from several ongoing surveys in 
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the Bootes field. These were used as training sets to adjust the templ ates and photonietric 



zero points to improve overall redshift accuracy and reliability (see iBrodwin et al.l 12006 

for details). Comparison with these spectroscopic samples shows that an rms dispersion of 

0"^ ~ 0.06(1 + z) is achieved for 95% of galaxies to at least z = 1.5. Subsequent follow- 



up spectroscopy of high redshift cand idate clusters ( §6.2t see also [Stanford et al.l 12005 






Brodwin et al.ll2006l : lElston et al.l 120061 ) verify this accuracy. 



A key output of the iBrodwin et al.l (120061 ) technique is the redshift probability function 
for each object, P{z), derived directly from the redshift-axis projection of the full redshift- 
SED likelihood surface. To transform these simple redshift likelihood functions into true 
probability distribution functions, a prior consisting of the observed redshift distribution was 
applied. This was measured in the FLAMEX region using high-quality BwRIJKs[3.6][4.5] 
photometric redshifts. Comparison to the spectroscopic sample illustrates that the resulting 
distribution functions are statistically valid in the sense that integrated areas accurately 
represent redshift probabilities at the 1, 2, and 3a levels. These P{z) functions are the input 
to the wavelet detection algorithm discussed in ^ and are also used in ^ 

Note that the P (z) distributions shows relatively little dependence on galaxy type 
( IBrodwin et al.l 120061 ) . The reason is that while all galaxies are selected to have at least 
5 sigma detections in [4.5], red galaxies are often only marginally detected or even unde- 
tected in our optical images, since they have very little blue light. This is particularly true 
at z > 1 where the fc-correction for early types is large in the optical. Blue galaxies have 
good i?^i?/[3.6][4.5] photometry, because they have more blue hght. Thus even though blue, 
late-type galaxy SED's have smaller breaks, they have more extensive (useful) photometry. 
These compensating effects lead to effectively type-independent photometric redshifts for 
z > 1 galaxies. 



5. Cluster Detection 



We employed a wavelet analysis to identify galaxy clusters within the Bootes region. 
Wavelet decomposition is a commonly used t e chnique for cluster ideii t ification in X-ray im- 



ages (for example , see IVikhlinin et al.l 1 19981 : IValtchanov et al.l l2004l : lAndreon et al.l 12005 



Renter et al.ll2005l ). where it provides an effective means of identifying extended sources in 
the presence of contaminating point sources. In principle, a similar analysis can be used 
with optical and infrared data sets, using galaxies rather than X-ray photons to identify 
extended sources. As is well known, galaxy number counts are more susceptible to projec- 
tion effects than is bremsstrahlung emission from the ICM. This issue is one that must be 
dealt with for all optical and infrared cluster searches, and consequently most such searches 
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for distant clusters make explicit assumptions abo ut the properties of the distant cluster 



population, such as a n assumed density profile (e.g., iPostman et al.lll996l:l01sen et al.lll999 



Scodeggio et al.lll999l) or the presence of a red sequence (JGladders &: Yee 



2000 



20051 ). The 



SDSS C4 technique (JMiller et al.ll2005l ) does not require a red sequence, but it does demand 
that the colors of cluster galaxies are similar to one another. 

A significant advantage of the Bootes data set is that the photometric redshifts permit 
such assumptions to be minimized. The full photometric redshift probability distributions 
P{z) were used to construct weighted galaxy density maps within overlapping redshift slices 
of width Az = 0.2, stepping through redshift space in increments of 6z = 0.1. For each 
galaxy the weight in the map corresponds to the probability that the galaxy lies within the 
given redshift slice. Weighting in this fashion de-emphasizes sources for which the redshift is 
poorly constrained. It is worth emphasizing that all galaxies were included in construction 
of the density map, and consequently cluster detection should be relatively independent of 
SED type and hence independent of morphology. Finally, cluster detection will be insensitive 
to the resolution of the density maps provided that the pixel size is small compared to the 
angular extent of cluster cores at all redshifts. A resolution of 12" (~ 100 kpc) per pixel 
was used, which satisfies this criterion while being sufficiently large to keep computational 
overhead manageable. 

Galaxy cluster candidates were detected within each redshift slice by convolving the 
density map with the wavelet kernel. We use a Gaussian difference kernel of the form 



k{r) 



1 



crt 



a^ 



where cxi = 400 kpc and a2 = 1600 kpc, and which crosses zero near r = 1 Mpc. The 
scale of the kernel is fixed in physical rather than angular units, preserving our ability to 
uniformly identify comparable systems at different redshifts. The precise physical values of 
a are subject to refinement, but the selected values effectively isolate overdensities on the 
scale of clusters or groups. 

Galaxy cluster candidates were detected in each redshift slice of these wavelet smoothed 
galaxy density maps using a simple peak-finding algorithm. To establish a consistent sig- 
nificance level for the candidates, 1000 bootstrap simulations were carried out within each 
redshift slice. The existing P(z) distributions, right ascensions, and declinations were re- 
peatedly shuffled, convolved with the wavelet kernel, and candidates detected to find the 
threshold corresponding to one false positive per redshift slice within the Bootes field. A list 
of detections above this significance threshold was generated for each redshift slice. Because 
the right ascensions and declinations were shuffled independently, the true correlated back- 
ground was not preserved. Consequently, the contamination rate may be somewhat higher 
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than the one false positive per redshift slice expected if we had preserved the correlated 
background. For the current analysis we accept the somewhat higher contamination rate 
in exchange for improved completeness, particularly for the highest redshift clusters. The 
majority of clusters are detected in multiple redshift slices (a natural consequence of sam- 
pling in slices separated by step sizes finer than the galaxy redshift uncertainties). Multiple 
detections with small separations in positions and redshift slices were considered to be a 
single cluster, which was assigned an estimated redshift and position corresponding to the 
slice with the highest statistical significance. A total of 335 candidates were found in redshift 
slices with centers from 0.1 to 1.9, including 98 in slices with z > 1. 

While the detection technique is on three dimensional overdensities, it is not immune to 
projection effects. The expected number of clusters in each cylindrical bin of length Az = 0.2 
and radius 1 Mpc (the radius of the wavelet detection kernel) was computed to quantify the 
extent of projection. This calculation included both the random expectation due to the 
observed nu mber density and the excess clusters expected due to the observed clustering of 



this sample (JBrodwin et al.ll2007l ). The projection rate is most significant for lower redshifts 
because of the large angular size corresponding to 1 Mpc, affecting ~ 10% of systems at 
z = 0.5 to ~ 20% at 2; = 0.1. The projection rate decreases to below ~ 4% at ^ = 1 and is 
negligible at 2; > 1.5. Overall we estimate that ~ 10% of the candidates may be spurious, 
including the contamination noted in the previous paragraph. 



6. Cluster Redshifts and Members 

The peak of the summed P{z) distribution at the cluster detection location was taken 
as the initial estimate of each cluster's redshift. To improve this estimate, individual objects 
within a 1 Mpc radius of the cluster which included the cluster redshift within the la range 
of their P{z) functions were considered candidate cluster members. A refined estimate (^est) 
of the mean cluster redshift was obtained from the peak of the summed P{z) distribution 
for these members. From this analysis, 104 candidates have Zest > 1- 

In Table [Hand Figures 2 - 13 we present 12 of the cluster candidates that were detected 
via the above criteria, and subsequently spectroscopically confirmed at the Keck Observatory 
to lie at ;z > 1 (see g63). Two of these, ISCS J1434. 1+3328 and ISCS J1429.2+3357, have 
Zest = 0.98 but the spectroscopic mean redshifts are z = 1.057 and 1.058 respectively. Hence 
we report a total of 106 clusters at 2; > 1, of which roughly 10% may be expected to arise by 
chance or from projection effects, for the reasons noted in ^ Column 1 of Table 1 provides 
the catalog number of each cluster. The catalog numbers increase with decreasing detection 
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significance ([|5])lj. Column 2 is the lAU designation for each cluster, based on the (J2000) 
coordinates of the detection given in columns 3 and 4. Column 5 provides the Zest value 
described above, and column 6 gives the mean redshift of the spectroscopically confirmed 
members. Column 7 provides the number of photometric redshift members of the cluster, 
defined as galaxies within 1 Mpc of the cluster center and with integrated P{z) > 0.3 in the 
range Zest ± 0.06(1 + -Zest)- These galaxies are used to calculate mean cluster colors in §7.21 
Column 8 reduces the number in column 7 by the number of galaxies which satisfy these 
criteria for each cluster redshift over the entire field, scaled to the 1 Mpc radius area for each 
cluster redshift. Column 9 gives the number of spectroscopically confirmed member galaxies 
in each cluster (see below). Column 10 gives the mean / — [3.6] color for the photometric 
redshift member galaxies, and column 11 provides the sum of their luminosities in [4.5] 
relative to an L* value which evolves according to the "red spike" model shown in Figure 
m corrected for the average over the field. Column 12 gives the luminosity in [4.5] of the 
brightest photometric redshift member galaxy relative to L*. 

We define a z > 1 cluster as spectroscopically confirmed if it contains at least 5 galaxies 
in the range ^est ± 0.06(1 + 2;ost) and within a radius of 2 Mpc, whose spectroscopic redshifts 
match to within ±2000(1 + {zgp)) km/s. The spectroscopic redshifts must also be of class 
A or B. Class A spectra have unambiguous redshift determinations, typically relying upon 
multiple well detected emission or absorption features. Class B spectral features are reliable 
but are less well detected. The radius threshold for spectroscopic membership is larger than 
the 1 Mpc used for photometric redshift members for practical reasons: most of the z > 1 
redshifts reported here were obtained with slitmasks extending out to approximately 2 Mpc 
from the cluster center. Typically < 10 photometric redshift members within 1 Mpc could 
be accommodated in the mask design, and often no redshift could be determined from the 
resulting spectra. 

The least significant cluster in Table 1 is ISCS J1434. 5+3427, which is the 327th most 
significant detection out of the sample of 335, and the 100th most significant detection at 
z > 1. Despite its relatively low detection significance, ISCS J1434. 5+3427 has a striking 
fjlar nentary morphology (Figure [9l), and has eleven spectroscopically confirmed members (see 



also iBrodwin et al.l l2006l ) 



■^These numbers extend beyond 335 for continuity with an earlier, preliminary version of the catalog used 
to plan spectroscopic and other foUowup observations. 



Table 1. z > 1 ISCS Clusters with Spectroscopic Confirmation^ 



Rank 


ID 


RA 


Dec 


^est 


\^sp/ 


Np 


Np(bc) 


Nsp 


a -[3.6]) 


J^totl J^ 


-^bcg/-^ 


152 


ISCS_J1434.1+3328 


14:34:10.37 


+33:28:18.3 


0.98 


1.057 


32 


13 


6 


4.83 


15 


2.5 


51 


ISCS_J1429.2+3357 


14:29:15.16 


+33:57:08.5 


0.98 


1.058 


45 


26 


7 


4.80 


34 


4.2 


19 


ISCS_J1433. 1+3334 


14:33:06.81 


+33:34:14.2 


1.02 


1.070 


57 


38 


20 


4.93 


45 


4.3 


123 


ISCS_J1433.2+3324 


14:33:16.01 


+33:24:37.4 


1.01 


1.096 


31 


11 


6 


4.94 


15 


4.6 


17 


ISCS_J1432.4+3332 


14:32:29.18 


+33:32:36.0 


1.08 


1.112 


49 


31 


23 


5.16 


47 


5.3 


34 


ISCS_J1426. 1+3403 


14:26:09.51 


+34:03:41.1 


1.08 


1.135 


31 


13 


7 


5.18 


26 


4.0 


14 


ISCS_J1426.5+3339 


14:26:30.42 


+33:39:33.2 


1.11 


1.161 


52 


35 


5 


5.15 


47 


2.8 


342 


ISCS_J1434.5+3427 


14:34:30.44 


+34:27:12.3 


1.20 


1.243 


27 


12 


11 


5.51 


27 


6.0 


30 


ISCS_J1429.3+3437 


14:29:18.51 


+34:37:25.8 


1.14 


1.258 


23 


7 


9 


5.30 


11 


6.5 


29 


ISCS_J1432.6+3436 


14:32:38.38 


+34:36:49.0 


1.24 


1.347 


30 


17 


8 


5.65 


31 


4.3 


25 


ISCS_J1434.7+3519 


14:34:46.33 


+35:19:33.5 


1.37 


1.373 


19 


9 


5 


5.77 


17 


3.9 


22 


ISCS_J1438.1+3414 


14:38:08.71 


+34:14:19.2 


1.33 


1.413 


25 


15 


10 


5.75 


24 


3.5 



^Columns are explained in 
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6.1. AGES Spectroscopy 

The extensive AGES spectroscopic database (C Kochanek et al. in preparation) can 
be used to spectroscopically confirm clusters at 2 < 0.5, a redshift at which L* corresponds 
to / = 20 (Figure 1), the magnitude hmit for AGES spectroscopy. Note that while AGES 
is 94% complete to / = 18.5, sparse sampling is used for fainter galaxies, and currently 
spectra are available for ~ 40% of galaxies with 18.5 </ < 20. Because of the high surface 
density of AGES redshifts and the larger angular scale aX z < 0.5, and the fact that the 
[4.5] fiux limit samples substantially less luminous galaxies at these redshifts, two additions 
to the criteria for spectroscopic confirmation used at 2; > 1 were imposed. To contribute to 
spectroscopic confirmation for cluster candidates with 2;est < 0.5, galaxies were required to 
be more luminous than L* + 1 in [4.5] for a red spike evolving model. Also, the average 
surface density of AGES galaxies over the Bootes field which would satisfy our high redshift 
confirmation criteria was calculated as a function of redshift, and 5 galaxies above this field 
level were required for a cluster to be confirmed. 

Of the 335 cluster candidates, 80 have Zest < 0.5. However seven lie near the edges 
of the Bootes field and hence were not well observed in AGES. Of the remaining 73, 61 
candidates are confirmed by AGES spectroscopy using the criteria just described. This 
criterion is perhaps overly stringent: it rejects two clusters at z = 0.2 with 17 or more 
matching AGES redshifts, but with only 4 (rather than 5) for galaxies brighter than L* + 1 
after field correction. All of the other 10 candidates observed by AGES which do not meet 
the confirmation criteria have z^st > 0.37, and half have Zest > 0.45. 

Given the sparse sampling of AGES at / > 18.5, which corresponds to L* at z > 0.3 
(Figure 1), we believe this confirmation rate validates our estimate that only ~ 10% of our 
cluster candidates arise by chance or due to projection effects. Further details on 2; < 1 
clusters will be provided in A. Gonzalez et al. (in preparation). 



6.2. Keck Spectroscopy 

Most of the high redshift spectroscopic confirmation of ISCS clusters has been obtained 
at Keck Observatory. Three clusters observed with Keck have been reported previously: 
ISCS J1432.4+3332 {z = 1. 112), ISCS J1434.5+3427 (z = 1.243), and ISCS J1438. 1+3414 



(z = 1 .413) are presented in lElston et al.l (120061 ) . iBrodwin et al.l (120061 ) . and [Stanford et al. 



( I2OO5I ). respectively. Here we provide spectroscopic confirmation for an additional nine 
clusters at z > 1, as well as some new spectroscopic information for the initial three clusters. 
Table [2] details new observations of ISCS clusters. Table [3] provides properties of previously 
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Fig. 2. — Composite NDWFS and Spitzer IRAC Bw, I, [4.5] color image of cluster 
ISCS J1434. 1+3328 at {zgp) = 1.057. North is up and east is left, and the field size is 
5' square (~ 2.4 — 2.5 Mpc for z = 1 — 1.5). The published version of these figures will label 
photometric redshift members, and spectroscopic members and non-members. 
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Fig. 3.— As for Figure El but for cluster ISCS J1429.2+3357 at (^^p) = 1.058. 
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Fig. 4.— As for Figure El but for cluster ISCS J1433.1+3334 at (^^p) = 1-070. 
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Fig. 5.— As for Figure El but for cluster ISCS J1433.2+3324 at (^.p) = 1.096. 
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Fig. 6.— As for Figure El but for cluster ISCS J1432.4+3332 at {z,p) = 1.112. 
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Fig. 7.— As for Figure El but for cluster ISCS J1426.1+3403 at (^.p) = 1.135. 
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Fig. 8.— As for Figure El but for cluster ISCS J1426.5+3339 at (^^p) = 1.161. 
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Fig. 9.— As for Figure El but for cluster ISCS J1434.5+3427 at (^^p) = 1.243. 
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Fig. 10.— As for Figure El but for cluster ISCS J1429.3+3437 at (^^p) = 1-258. 
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Fig. 11.— As for Figure El but for cluster ISCS J1432.6+3436 at {z,p) = 1.347. 
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Fig. 12.— As for Figure El but for cluster ISCS J1434.7+3519 at (^^p) = 1-373. 
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Fig. 13.— As for Figure El but for cluster ISCS J1438. 1+3414 at {z,p) = 1.413. 
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unreported, spectroscopically confirmed cluster members, in decUnation order for IRAC 
sources, followed by serendipitous sources. Figure [H] shows three example spectra obtained 
in April 2007 with Keck/DEIMOS. 



6.2.1. LRIS Observations 
We obtained deep optical slitmask spectroscop y for several clu sters using the dual- 



beam Low Resolution Imaging Spectrograph (LRIS; lOke et al.lll995l ) on the 10 m Keck I 
telescope during three observing runs between 2005 and 2007. Slitmasks generally included 
approximately 15 objects with photometric redshifts consistent with cluster membership and 
within 4 arcmin of the nominal cluster centers. Additional IRAC 4.5/im selected sources were 
included to fill out the slitmasks. Slitlets had widths of 1.3 arcsec and minimum lengths of 10 
arcsec. We employed the D580 dichroic which splits the light at ~ 5800 A between the two 
channels of LRIS. On the red side, the 400 fine grating, blazed at 8500 A, was used to cover 
a nominal wavelength range of 5800 to 9800 A, varying somewhat depending on the position 
of a slit in the mask. For objects filling a slitlet, the spectral resolution for this instrument 
configuration is ~ 9 A (i? ~ 900), as determined from arc lamp spectra. On the blue side, 
the 400 line grism, blazed at 3400 A, provided coverage from the atmospheric cutoff (~ 3200 
A) up to the dichroic cut off. For objects filling a slitlet, the spectral resolution for this 
instrument configuration is ~ 8 A (i? ~ 450). 

We obtained multiple exposures for each mask, usually with 1800 s per individual ex- 
posure. Table [2] details the exposure times and observing conditions. The observations were 
carried out with the slitlets aligned close to the parallactic angle, and objects were shifted 
along the long axis of the slitlets between exposures to enable better sky subtraction and 
fringe correction. 

The slitmask data were separated into individual spectra and then reduced using stan- 
dard longslit techniques. The multiple exposures for each slitlet were reduced separately and 
then coadded. The spectra were reduced both without and with a fringe correction; the for- 
mer tends to yield higher quality object spectra at the shorter wavelengths, while the latter 
is necessary at the longer wavelengths. Calibrations were obtained from arc lamp exposures 
taken immediately after the object exposures for the red side, and from arc lamp exposures 
taken during the afternoon for the blue side. Corrections for small offsets in the wavelength 
calibration were obtained by inspection of the positio ns of sky lines in the objec t spectra. 



Using longslit observations of the standard stars from iMassey fc Gronwalll (Il990l ) obtained 
during the same observing runs, we achieved relative flux calibration of the spectroscopy. 
While slit losses for resolved sources preclude absolute spectrophotometry from the slit mask 



Table 2. Keck and Subaru Observations of z > 1 ISCS Galaxy Clusters 



Cluster 



(^sp) 



Inst. 



UT Date 



Exp. Time (s) 



Conditions 



ISCS_J1434.1+3328 
ISCS_J1429.2+3357 
ISCS_J1433. 1+3334 
ISCS_J1433.2+3324 
ISCS_J1432.4+3332 



ISCS_J1426. 1+3403 
ISCS_J1426.5+3339 

ISCS_J1434.5+3427 

ISCS_J1429.3+3437 
ISCS_J1432.6+3436 
ISCS_J1434.7+3519 
ISCS_J1438. 1+3414 



1.057 


DEIMOS 


2007 Apr 18-19 


6x1800 


clear; 0'.'8 - 2'.'0 


1.058 


DEIMOS 


2006 Apr 26 


4x1500 


not photometric 


1.070 


DEIMOS 


2007 Apr 19 


3x1200 


clear; 0'.'8 


1.096 


DEIMOS 


2007 Apr 18-19 


6x1800 


clear; ^'8 - 2'.'0 


1.112 


LRIS 


2005 Jun 03 


5x1800 


clear; 0'.'9 (Elston et al. 2006) 




FOCAS 


2006 Apr 21 


5x1200 






DEIMOS 


2007 Apr 19 


3x1200 


clear; 0'.'8 


1.135 


LRIS 


2007 May 19 


3x1800 


clear; 0'.'9 


1.161 


LRIS 


2006 Apr 04 


4x1800 


clear 


1.243 


FOCAS 


2006 Apr 22 


5x1200 






DEIMOS 


2007 Apr 18 


8x1800 


likely cirrus; I'.'O - 1'.'5 


1.258 


LRIS 


2006 Apr 05 


3x1200 




1.347 


LRIS 


2007 May 21 


7x1800 


mostly clear; 0'.'8 


1.373 


LRIS 


2005 Jun 02 


7x1800 




1.413 


FOCAS 


2006 Jun 28 


9x1200 






DEIMOS 


2007 Apr 19 


7x1800 


clear; 0'.'9 



to 

CO 
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Fig. 14. — Optical spectra of three z > 1 members of IRAC Shallow Survey clusters. The 
spectra were obtained with Keck/DEIMOS in April 2007 and flux calibrated using an archival 
sensitivity function. All three sources have class A spectra (§6]) and show clear Ca H and 
K absorption, typical of early-type galaxies. The bottom source also shows strong [O II] 
emission, indicative of star formation. All three sources have I ~ 21.7. Imperfect sky 
subtraction is evident in the top spectrum near 7300 A, and the 7600 A atmospheric A-band 
is marked in all 3 panels. 
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data, the relative calibration of the spectral shapes should be accurate. One-dimensional 
spectra were extracted from the sum of all the reduced data for each slitlet for both the red 
and blue sides. For the targets in the high-redshift clusters, generally only the red side data 
proved useful. 



6.2.2. DEIMOS Observations 
Additional sp ectroscopy was ob tained with the Deep Imaging Multi-Object Spectro- 



graph (DEIMOS; iFaber et al.l 120031 ) on the 10 m Keck II telescope, a second generation 
instrument with significantly more multiplexing capabilities as compared to LRIS, albeit 
without the blue sensitivity. During an observing run in April 2006, we targeted Boote s ac- 



tive galaxy candidates selected on the basis of mid-infrared colors (e.g.. lStern et al.ll2005l ). but 



included candidate cluster members in one mask. In April 2007, while observing host galax- 
ies of high redshift cluster supernovae (Perlmutter et al., in preparation), we also observed 
candidate cluster members, typically with more than one high redshift cluster candidate 
observed in each wide area mask. 

For both observing runs, the 600ZD grating (Abiazo = 7500 A; AAfwhm = 3.7 A) and a 
GG455 order-blocking filter were used. DEIMOS data were processed using a slightly modi- 
fied version of the pipeline developed by the DEEP2 team at UC-BerkelejQ. Although neither 
run was completely photometric, relative flux c alibration was achieved from observations of 



standard stars from lMassey &: Gronwalll (Il990l ). 



6.3. Subaru Spectroscopy 

FoUowup spectrosc opic observations were also obtained with the FOCAS spectrograph 



(iKashikawa et al.l |2002| ) on the 8.2-m Subaru telescope in April and June 2006. For these 
observations, we used the instrument in multi-object spectroscopy mode with 0.8 arcsec 
width slits, 300R grism, and S058 order-sorting filter, which provided about 15 spectra 
from 5800A to lOOOOA with a spectral resolution R ~ 500. Total exposure times were 2-6 
hours. The FOCAS data were reduced with IRAF using standard methods including a fringe 
correction. Wavelength calibration was done using OH airglow emission lines. Absolute flux 
was calibrated using standard star (Feige 34, Wolf 1346, Hz 44) spectra taken during the 
same nights. 



http : //astro . berkeley . edu/$\sim$cooper/deep/spec2d/ 
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6.4. Notes on Individual Clusters 

6.4.1. ISCS JI434. 1+3328 (z = 1.057) 

ISCS J1434. 1+3328 was observed with DEIMOS on a mask optimized for observing a 
high redshift supernova identified in a different ISCS high redshift cluster (Perlmutter et 
al., in preparation) which has not yet been spectroscopically confirmed. Four candidate 
members of ISCS J1434. 1+3328 were targeted, and two more spectroscopic members were 
identified from additional slitlets targeting 4.5 /xm selected sources to fill the mask. All six 
confirmed members have clearly identified Ca H and K absorption lines and D4000 breaks, 
and none of the galaxies exhibit emission features in the wavelengths covered by the DEIMOS 
observations {^ 5500 — 9500 A). 

6.4.2. ISCS JI429. 2+3357 (z = 1.058) 

ISCS J1429. 2+3357 was observed with DEI MOS, on a slitmas k optimized for observing 



faint, mid-infrared selected AGN candidates (e.g.. lStern et al.ll2005l ). Eight candidate cluster 
members were included in the slitmask, six of which were confirmed spectroscopically to re- 
side at z^ 1.06. Only the brightest galaxy, IRAC J142912.9+335808, shows [O II] emission; 
the remaining redshifts were derived on the basis of Ca H and K absorption lines and/or 
D4000 breaks. In addition, one galaxy targeted spectroscopically as an IRAC-selected AGN 
candidate, IRAC J 1429 16. 1+335537, is a cluster member, bringing the tally to seven spec- 
troscopically confirmed cluster members. The spectrum of this source shows strong, narrow 
(400 km s~^) emission lines from [Ne V] A3426, indicating that it is indeed an active galaxy. 



6.4.3. ISCS JI433. 1+3334 (z = 1-070) 

ISCS J1433. 1+3334 was observed with DEIMOS on a mask optimized for observing a 
high redshift supernova identified in ISCS J1432. 4+3332 (Perlmutter et al., in preparation). 
Eight candidate members of ISCS J1434. 1+3328 were targeted, of which six were confirmed 
as cluster members, one was found to be slightly foreground to the cluster, and one was 
slightly behind the cluster. All eight candidates show clear D4000 breaks. Many additional 
cluster members were identified on this mask, either serendipitously or as targeted IRAC- 
selected, z > 1 galaxies, bringing the total number of spectroscopically confirmed cluster 
members to 20. As seen in the bottom two panels of Figure [TH the confirmed sources show 
a range of spectral properties. While most show Ca H and K absorption and D4000 breaks. 



-33- 



some also show emission features likely due to either star formation or AGN activity. 



6.4.4. ISCS JI433. 2+3324 (z = 1-096) 

ISCS J1433.2+3324 was observed with DEIMOS on the same mask as ISCS J1434.1+3328. 
Five candidate members of ISCS J1434. 1+3328 were targeted, of which two were confirmed 
as cluster members, one was found to be foreground, and two yielded inconclusive spec- 
tra. Four additional cluster members were identified in the same mask, from IRAC-selected 
sources. The two confirmed cluster members which were specifically targeted show strong 
Ca H and K absorption and lack emission lines. The other four confirmed members all show 
[O II] emission, with three of the four also showing D4000 breaks. 



6.4.5. ISCS JI432. 4+3332 (z = 1.112) 

The sp ectroscopic observa tions which confirmed ISCS J1432. 4+3332 at z = 1.11 are 
d escribed in Elston et al.l (120061 ). but no data on individual sources was presented there. The 
Iston et al.l (120061 ) result was based on nine spectroscopically confirmed cluster members, 
two of which were selected not as cluster members, but rather as mid-IR selected AGN, and 
that data is now included in Table [3] of this paper. One of the candidate members of this 
cluster hosted a high-redshift supernova (Perlmutter et al., in preparation) and was thus 
targeted for additional DEIMOS and FOCAS slitmask spectroscopy. An example DEIMOS 
spectrum for this cluster is shown in the top panel of Figure [HJ In total, there are now 23 
spectroscopically confirmed members in the cluster. 



6.4.6. ISCS JI426. 1+3403 (z = 1.135) 

Seven candidate members of this cluster were confirmed spectroscopically during our 
LRIS observations in May 2007. All seven galaxies show red continuum emission and/or a 
clear D4000 break. Two of the sources also show [O II] emission. 



6.4.7. ISCS JI426. 5+3339 (z = 1.161) 

Four of the five spectroscopically confirmed members of this cluster show [O II] emission, 
a higher proportion than is typical for this program. This is possibly a selection effect; sources 



-34- 



with line emission are the easiest to spectroscopically confirm. All five confirmed members 
show breaks typical of early- type galaxies (e.g., D2900 and/or D4000). 



6.4.8. ISCS JI434. 5+3427 (z = 1.243) 



This cluster is discussed in iBrodwin et al.l ( l2006l ) , where eight spectroscopic members 



were presented. Table [3] of the current paper does not duplicate those data, and instead lists 
three additional spectroscopic members that have since been identified. 



6.4.9. ISCS JI429. 3+3437 (z = 1.258) 

This cluster has nine spectroscopically confirmed members. One is an optically-bright 
AGN from the AGES survey, while the rest were confirmed spectroscopically by LRIS. Four 
of the Keck/LRIS sample show [O II] emission, while the other four show only spectral breaks 
and absorption lines. 



6.4.10. ISCS J 1432. 6+34 36 (z = 1.34 7) 

This cluster has eight spectroscopically confirmed members, all from our Keck/LRIS 
observations in May 2007. Only one of the sources shows [O II] emission; the rest of the 
redshifts are on the basis of continuum breaks at 2640 and 2900 A, as well as Mg II A 2800 
absorption. 



6.4.11. ISCS JI434. 7+3519 (z = 1.373) 

This cluster has five confirmed members, one of which was serendipitously identified in 
the spectroscopy, all from our Keck/LRIS spectroscopy in June 2005. Three of the members 
show [O II], and the other two show spectral breaks characteristic of early- type galaxy 
spectra. 
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6.4-12. ISCS J1438. 1+3414 (z = 1.413) 



This cluster was first published in IStanford et al.l (120051 ) . at which time it was the 
highest redshift galaxy cluster known. Two candidate cluster members hosted supernovae 
in the HST/ ACS program of Perlmutter et al. (in preparation), and this field thus has been 
the target of additional spectroscopy from Subaru and Keck. Five new cluster members have 
been confirmed , listed in Table [31 Data on the original five members is not duplicated from 
Stanford et al.l (120051 ). Of the five new members, one is an [O II] emitter, serendipitously 
identified in a slitlet targeting another source. Three show only Ca H and K absorption 
lines, with no emission lines identified. IRAC J143816. 8+341440 (22.3 in Table 3) is an 
AGN, showing [Ne IV] , [O II] , and [Ne III] emission lines. 



6.5. Cluster Masses 



With 20 or more spectroscopic redshifts, it is possible to estimate cluster masses via 
scaling relations using the velocity dispersion. The line of sight velocity dispersion for 
the 20 spectroscopic member galaxies in cluster 19 (ISCS J1433. 1+3334) at z = 1.070 is 
760 km s~^ in the rest frame, and for the 23 spectroscopic member galaxies in cluster 17 
(ISCS J1432.4+3332) at ;z = 1.112 it is 734 km s'^ 



of' 



For these two clusters, which are among the richest in the sample, the velocity dispersion 
■ 750 km s~^ correspond s to a virial mass of ~ 3.8 r^ x IO^^Mq where r^ is the virial radius 

T he x-r ay temperature corresponding to 



in Mpc (e.g. equation 4 of ICarlberg et a. 
a = 750 km s"^ is 3.9 keV 



4 X lO^^Mo (Shimizu et al. 



able 4 of 



19961 



Xue fc Wul I2OOOD , which in turn gives a mass of 



20031) . This is consistent with the average halo mass of 



IO^^Mq estimated by iBrodwin et al.l (120071 ) for the sample. 



Stellar luminosities can also be used to make a rough estimate of cluster masses, scaling 
to the Com a cluster via the r ed sp ike model. Within the 650 x 850 kpc region of Coma 
sampled by lEisenhardt et al.l (120071 ) the integrated K-hand luminosity to L* + 1 is 54L*, 
slightly higher than the values shown in column 11 of Table [H but in a significantly smaller 
effective radius of 0.42 Mpc. In their Figure 1, Geller, Diaferio, & Kurtz (1999) show a mass 
for the Coma cluster of ~ 3.5 x IO^^Mq within this radius (for h=0.7). At r=l Mpc they 
show about double that mass, and suggest the total mass for Coma is about double this 
again. Hence if the profile and M/Lk for our clusters is similar to Coma, allowing for red 
spike model evolution in L^, the Ltot/L* values in 1 Mpc radius provided in Table 1 scale 
to total cluster masses of ~ 1 — 6 x IO^^Mq (i.e. ~ 0.1 — 0.4Mcoma)- Clusters 17 and 19 are 
at the top of this range, again providing mass estimates consistent with those found using 
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the velocity dispersion. More detailed exploration of this approach will require photometry 
which takes account of the extent to which a 5" diameter aperture fails to include all of the 
light from such galaxies, or includes light from multiple galaxies (as determined from higher 
spatial resolution imaging). 



Discussion 



As noted in the introduction, out to 2; ~ 1 substantial evidence exists which is consistent 
with an extremely simple formation history for cluster galaxies, in which their stars are 
formed in a short burst at high redshift, and they evolve quiescently thereafter (ie a "red 
spike" model). The colors of luminous cluster galaxies out to ^ ~ 1 typically fall on a tight 
sequence which is red relative to field galaxies at similar redshift (the "red sequence"), with 
a mean color which evolves as the red spike model predicts. 

The most luminous red sequence galaxies are the reddest, a correlation whi ch is at- 
trib uted to higher metalli city in higher mass galaxies (JKodama fc Arimotd 119971 : but see 
also iFerreras fc Silkl |2003| ). This correlation (i.e the color-magnitude or mass- metallicity 
relation) is explained by additional cycles of star formation and enrichment in more mas- 
sive ga laxies, as they are able t o retain their gas more effectively against super novae- driven 
winds ( lArimoto fc Yoshiil 119871 ). In the red spike formation paradigm, this is a natural con- 
sequence. In the context of the hierarchical merging galaxy formation models, a unique 
correlation between stellar mass and metallicity is a greater challenge. While the inclusion 
of feedback (whet her by supernovae or by AGN) in hierarchical models s tops the buildup of 
galaxy mass (e.g. iKauffmann fc Charlotlll998l : iNagashima fc Yoshiil |2004J ). such models have 
difficulty reproducing the exact color and slope of the color-magnitude relation in clusters 
over the full redshift range for which it has been measured. Recent work shows promise, 
how ever, as feedback beha vior and other "gastrophysical" effects are taken more into account 
(e.g. 



De Lucia et al.ll2006f ). 



If the onset of the star formation "spike" is simultaneous for all cluster galaxies, and if the 
color-magnitude relation is caused by more protracted spikes in more massive galaxies, then 
as one approaches the star-forming epoch, massive galaxies might be expected to eventually 
become bluer, reversing the slope of the color-magnitude r elation. In fact no measurable 



chang e is seen in color-magnitude slope out to 2; ~ 1 (e.g [Stanford et al.l Il998l : iMei et al. 
20061 ). which suggests either formation redshifts well before z = 1, or that the spikes begin 
earlier in more massive galaxies, perhaps ending rather than beginning simul taneously. The 



small and unchanging scat t er of the red sequ ence in clusters out to z ~ 1 (e.g. [Stanford et al. 



19981 : iBlakeslee et al.ll2003l : iTran et al.ll2007l ) also argues for synchronized or very early spikes 
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and against a primarily age-based origin of the color- magnitude relation, since the scatter 
would increase by a much larger amount, and faster, than has been observed. 

Testing to what extent red spike models remain consistent with cluster galaxy data at 
z > 1 is one of the primary motivations for the present study. We begin by constructing the 
color-magnitude relations for the 12 spectroscopically confirmed clusters in Table [H 



7.1. Color-Magnitude Diagrams 

Figures [15] and [12] present color magnitude diagrams for the clusters listed in Table 
[1] The / — [3.6] color was selected because these filters bracket the 4000A break most 
tightly at 2; > 1. The symbol area is proportional to the integral of the object's redshift 
probability distribution over the range z^st =1= 0.06(1 -1- Zest), which is the rms dispersion in 
individual photometric redshifts ( Brodwin et al.ll2006l ). Circled symbols are spectroscopically 



confirmed members, while crosses indicate objects known not to be members on the basis 
of spectroscopy. Figures [15] and [16] include objects within 1 Mpc of the cluster center, with 
> 5cr detections in both [3.6] and [4.5], and with > 2a detections in /; or with spectroscopic 
redshifts. As noted in §3.2[ the 5a limit in [3.6] is 18.6 mag, so the limiting factor in object 
selection for Figures [T5] and [T6] is from the 5a limit of 17.8 mag in [4.5]. We have used the 
0.1 Gyr burst, Zf = 3 red spike model shown in Figure [T] to calculate the corresponding [3.6] 
mag as a function of redshift, and this is shown in Figures [T5] and [T6] by the vertical dotted 
lines. The I limit is shown by the diagonal dotted lines in Figures [15] and [16] The vertical 
dashed line plots the expected L* magnitude at [3.6] for the red spike model. The sloped 
dashed hne sho ws the observed U — ij" color- magnitude slope of 0.22 for the Coma cluster 



from the data of lEisenhardt et al.l (120071 ). normalized at the red spike model [3.6] magnitude 
and / — [3.6] color for an L* galaxy. In this redshift range observed I — [3.6] is close to 
rest U — H. Note that while the normalization evolves according to the red spike model, a 
constant, unevolving slope value is used. 

The brightest galaxy likely to be a member is typically 1 to 2 magnitudes brighter 
than the expected L* magnitude for the red spike model (see Table H]). It is noteworthy 
that the brightest galaxy in the Coma cluster is ~ 6 times brighter than L* in the H band 



( Ide Propris et al.l ll998l : lEisenhardt et al.l 120071 ) , suggesting less than a factor of two growth 



in stellar mass in such galaxies since z ~ 1.5. HST imaging should be used to assess this 
more carefully. 

Because these clusters were not selected on the basis of containing a red sequence (^, 
the fact that the highest probability cluster members tend to track the passively evolving 



-38- 



CO 



7 
6.5 

6 
5.5 

5 
4.5 

4 
3.5 

3 
2.5 

7 
6.5 

6 
5.5 

5 
4.5 

4 
3.5 

3 
2.5 

7 
6.5 

6 
5.5 

5 
4.5 

4 
3.5 

3 
2.5 



\ \ I p" 

152(z=1.057) I 



® 






* ®~!S)f4 



*« 



X X X ' v'^ 



^ 



^i^Hh^'^ 



t»i 



I I 1 I ■-!.. I 

123(z=1.096') I ''■{-. 

' lllllliJ 

• ■■■ « rr li 




15.5 16 16.5 17 17.5 18 18.5 19 15.5 16 16.5 17 17.5 18 18.5 19 



[3.6] 



■ /P{z) = 0.3 
O Spectroscopic member 
X Spectroscopic nonmember 



Fig. 15. — Observed / — [3.6] vs. [3.6] for galaxies within 1 Mpc of the center of the 
spectroscopically confirmed clusters at z > 1, in order of increasing redshift. The catalog 
number and spectroscopic redshift for each cluster from Table [1] are given in the upper left 
corner of each sub-panel. The symbol area is proportional to the galaxy's integrated redshift 
probability in the range z^gt ± 0.06(1 + -Zest)- A symbol with an integrated probability of 0.3 
is shown for reference in the lower right portion of the figure. Additional details are provided 
ingm 
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Fig. 16. — Observed / — [3.6] vs. [3.6] for the remainder of the spectroscopically confirmed 
clusters at 2; > 1, in order of increasing redshift. See figure [15] for details. 
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Coma cluster sequence shown by the dashed hne is significant. Evidently the red sequence 
persists in dense environments to z = 1.4, even when not used as a selection criterion. In 
addition there is no indication that the slope of the color-magnitude relation has changed 
sign, and in fact a non-evolving Coma cluster slope appears consistent with the data. 

Comparing the data for e.g. cluster 14 (ISCS J1426.5+3339 at z = 1.16) vs. 29 (ISCS 
J1432. 6+3436 at z = 1.34) suggests real differences in the scatter of the color-magnitude 
relation do exist. The presence of luminous galaxies substantially bluer than the red sequence 
in several clusters, a number of them spectroscopically confirmed, is also noteworthy. These 
may represent a population of massive star-forming gal axies which fade onto the red sequence 



in rich clusters, but is not found in field surveys (e.g. iBell et al.l 120041). It shou l d be n oted 



that the very lo w scatters reported for z > 1 clusters in e.g. iBlakeslee et al.l (120031 ) and 



Mei et al.l (120061 ) are calculated for galaxies known to have early-type morphologies, using 
much deeper photometry than the survey/ discovery data used here. Further investigation 
of the color-magnitude relations is deferred to a future paper which will use deeper Spitzer 
and HST imaging presently being obtained. 

Although the survey data used to identify the clusters does not enable accurate pho- 
tometry of individual galaxies in clusters, we can calculate mean properties for galaxies in 
each cluster with some confidence. In the next section we briefiy explore the color evolution 
and color-magnitude relation of the entire cluster sample. 



7.2. Color vs. Redshift 

Figure [T7I plots the average J — [3.6] color for galaxies within 1 Mpc of the cluster centers, 
and whose integrated redshift probability distribution in the range 2;est±0.06(l + Zest) exceeds 
0.3. The mean values are calculated after iteratively clipping 3cr outliers. As implied by the 
color-magnitude diagrams for the spectroscopically confirmed clusters ( §7.11) . the red spike 
model (0.1 Gyr burst at Zf = 3) provides a remarkably good fit to the empirical data, 
validating the survey design assumptions illustrated in Figure [TJ 

Much of the increase in / — [3.6] with redshift is due to the change in the rest frame 
wavelengths observed, and in Figure [18] we plot the average color relative to a no-evolution 
model (i.e. the k-corrected color). The no-evolution model is simply the red spike model 
at z=0, i.e. at an age of 11.3 Gyr for the stellar population. Clearly the evolving red spike 
model is a better fit to the data than is the k-correction alone. Note that the photometric 
redshifts (§!]) are calculated using non-evolving templates, so this result is not a foregone 
conclusion, particularly at z > 1. 
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Fig. 17. — Mean observed / — [3.6] vs estimated cluster redshift Zest for galaxies within 
1 Mpc of the cluster center and with integrated redshift probabilities > 0.3 in the range 
Zest ±0.06(1 + ^est)- The sohd red line is the predicted /— [3.6] color for the red spike model. 
Filled symbols show the mean I — [3.6] colors for the spectroscopically confirmed z > 1 
clusters listed in Table [H 
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Fig. 18. — As for Figure [T71 but with the k-correction subtracted. Red squares are the mean 
observed / — [3.6] for galaxies brighter than the expected [3.6] magnitude for an L* galaxy 
for the red spike model. The solid red line is the predicted I — [3.6] color for the red spike 
model. Filled symbols show the mean I — [3.6] colors for the spectroscopically confirmed 
z > 1 clusters listed in Table [D 
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The mean /— [3.6] colors in Figure [T7| include galaxies down to the 13.3yuJy survey limit 
at 4.5/im, and therefore less luminous galaxies contribute to the mean at lower redshifts. The 
color-magnitude relation shows that less luminous galaxies are bluer, and massive galaxies 
are quite quiescent today, so a constant flux limit will lead to a systematic bias towards more 
luminous, massive galaxies and redder mean color with increasing redshift. On the other 
hand the excellent fit of the red spike model implies that a constant luminosity limit for the 
sample would select galaxies with smaller stellar masses at high redshift, since a constant 
stellar mass, passively evolving, will be more luminous as one approaches the formation 
redshift. To avoid this bias, in Figure [18] we also show (with red squares) the mean colors 
for galaxies brighter than the red spike passively-evolving L* in [3.6]. 

It is evident that the mean color of the more luminous galaxies is systematically redder, 
and hence that the color-magnitude relation continues to hold out to 2; ~ 1.5. The color 
offset is 0.1 magnitudes, independent of redshift. The persistence of the essentially unchanged 
color-magnitude relation slope, with a passively evolving intercept, out to lookback times 
within 4 Gyr of the Big Bang is a phenomenon that models of cluster galaxy formation 
must account for. The implication is that the correlation between high stellar population 
metallicity and stellar mass is already in place at z ~ 1.5, and that the star formation era 
remains well in the past for these galaxies. 

In Figure [19] the mean / — [3.6] colors are plotted with the predicted / — [3.6] color for 
the red spike model subtracted. The mean colors are for galaxies more luminous in [3.6] than 
the red spike L*. Note the red spike model is a good match to L* galaxy colors in figures [T7] 
and [IHl so the offset in color out to 2; ~ 1 is attributable to the color-magnitude relation. 

Other effects may be in play. In Table 1 the z^st values tend to be lower than the 
(zsp) values, which may lead to colors appearing ~ 0.1 mag redder relative to the red spike 
model than if (z^p) were available for all 2; > 1 clusters. The average colors of cluster 
members may still be redder than shown, however, as no attempt has been made to correct 
the average cluster colors for field contamination, and field galaxies tend to be bluer than 
cluster galaxies. Blending issues could lead to colors which are systematically biased. Given 
these uncertainties, the excellent agreement with the red spike model is the more remarkable. 
Nevertheless, at 2; > 1 the trend is towards colors which are increasingly redder than the 
red spike model, for both the full sample and the spectroscopically confirmed subset, though 
there is clearly a range. Hence even higher formation redshifts (as high as Zf = 30) are 
favored for most z > 1 clusters, within the context of red spike models. 
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Fig. 19. — As for Figure [TTl but with predicted I — [3.6] color for the red spike model 
subtracted. Means are shown only for galaxies brighter than the expected [3.6] magnitude 
for an L* galaxy for the red spike model. The black horizontal line is the 0.1 Gyr burst, 
Zf = 3 (red spike) model prediction. The purple, blue, green, orange, red, and black lines 
show the expected relative I — [3.6] color for a 0.1 Gyr burst starting at Zf = 1.5, 2, 4, 5, 30, 
and for no evolution (k-correction) respectively. Representative values of the uncertainty in 
each cluster's mean / — [3.6] color are shown at several redshifts. Filled symbols show the 
mean / — [3.6] colors for the spectroscopically confirmed z > 1 clusters listed in Tabled! 
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7.3. Clusters at 1.5 < z <2 

Figure [T71 shows that that there are relatively few cluster candidates in the IRAC sample 
with Zest > 1-5. Why is this? The red spike model fits, circumstantial evidence from 
BCG luminosity, and the persistence of the color-magnitude relation suggest that the stellar 
populations were formed - and assembled - well before z ~ 1.5. From Figure [U such clusters 
should be detectable by the IRAC Shallow Survey. 

A potential selection effect against 2; > 1.5 galaxies is that photometric redshifts in this 
range have increasingly broad redshift probability distributions, due in part to insufficient 
depth in photometry at shorter wavelengths. Simulations indicate a factor of two reduction in 
IRAC photometric error can help compensate, leading to a similar reduction in photometric 
redshift error in this redshift range. With the Spitzer Deep Wide-Field Survey legacy program 
now underway (PI D. Stern), such data will be available in the near future. The tighter 
redshift probabihty distributions would improve the contrast of high redshift clusters over 
the field, allowing them to meet our detection threshold. 

But the lack of massive z > 1.5 cluster detections may simply refiect a real decline in the 
space density of such objects. Models for the hierarchical growth of structure predict only 
about 1 cluster with z > 1.4 above 10^'^Mq for every 5 clusters with z > 1. Full exploration 
of this possibility will require more careful simulation and assessment of our observational 
selection effects. 



8. Summary 

We have identified 335 galaxy cluster and group candidates from a 4.5/im selected 
sample of galaxies in the IRAC Shallow Survey. Candidates were identified by searching 
for overdensities in photometric redshift slices, and 106 clusters are at 2; > 1. Roughly 10% 
of these candidates may be expected to arise by chance or from projection effects. To date, 12 
clusters have been spectroscopically confirmed at z > 1, as have 61 of the 73 clusters observed 
with AGES at 2; < 0.5. For the two z > 1 clusters with 20 or more spectroscopic members, 
total cluster masses of several lO^^M© are indicated, and the total mass estimated from the 
stellar luminosity yields comparable values. Color-magnitude diagrams in J — [3.6] vs. [3.6] 
for the z > 1 spectroscopically confirmed clusters reveal that a red sequence is generally 
present, even though clusters were not selected for this. The brightest probable member 
galaxy (at the spatial resolution of IRAC) in the spectroscopically confirmed z > 1 clusters 
remains 1-2 mag brighter than the passively-evolving L* luminosity. For the full cluster 
sample, the mean color of brighter galaxies within each cluster is systematically redder than 
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the mean color of all probable cluster member galaxies, implying that the mass-metallicity 
relation is already in place at z ~ 1.5. The mean / — [3.6] color of probable cluster members 
is well fit by a simple model in which stars form in a 0.1 Gyr burst beginning at Zj = 3. 
At 2; > 1 there is a tendency for mean cluster colors to favor formation redshifts Zf > 3, 
although a few are consistent with Zf ~ 2. This adds to the large body of evidence that 
galaxies in clusters were established at extremely early times. 
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Table 3. z > 1 Spectroscopic Cluster Members 



ID^ 


4.5]^ 


phot-z 


spec-z 


ISCS_J1434.1+3328 (^sp) 


= 1.057 


152.1 


15.45 


1.01 


1.057 


152.2 


15.99 


0.91 


1.055 


152.3 


17.00 


0.96 


1.054 


152.4 


17.47 


1.04 


1.065 


152.5 


16.40 


0.99 


1.055 


152.6 


16.76 


0.96 


1.055 


ISCS_J1429.2+3357 (^sp) 


= 1.058 


51.1 


15.10 


0.98 


1.059 


51.2 


17.06 


0.95 


1.056 


51.3 


17.83 


1.04 


1.055 


51.4 


17.16 


1.00 


1.060 


51.5 


17.01 


1.03 


1.059 


51.6 


17.19 


1.03 


1.054 


51.7 


16.08 


1.01 


1.060 


ISCS_J1433.1+3334 (^sp) 


= 1.070 


19.1 


17.49 


0.99 


1.075 


19.2 


16.67 


0.99 


1.075 


19.3 


17.04 


1.06 


1.076 


19.4 


16.00 


1.06 


1.066 


19.5 


16.89 


1.09 


1.067 


19.6 


16.47 


0.94 


1.065 


19.7 


17.70 


1.00 


1.063 


19.8 


17.15 


1.10 


1.079 


19.9 


16.75 


0.98 


1.074 


19.10 


16.30 


1.07 


1.066 


19.11 


16.21 


1.03 


1.066 


19.12 


16.66 


1.03 


1.075 


19.13 


15.44 


1.06 


1.064 


19.14 


17.03 


1.28 


1.076 


19.15 


17.51 


1.00 


1.065 
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Table 3 


— Continued 


ID^ 


4.5]^ 


phot-z 


spec-z 


19.16 


17.01 


1.27 


1.069 


19.17 


17.36 


1.04 


1.066 


19.18 


16.27 


0.86 


1.063 


19.19 


16.68 


2.45 


1.074 


19.20" 


- 


- 


1.067 


ISCS_J1433.2+3324 (^sp) 


= 1.096 


123.1 


15.32 


1.07 


1.094 


123.2 


16.45 


1.09 


1.093 


123.3 


16.68 


1.10 


1.094 


123.4 


16.99 


1.11 


1.107 


123.5 


16.64 


0.93 


1.094 


123.6 


17.35 


0.99 


1.091 


ISCS_J1432.4+3332 (^sp) 


= 1.112 


17.1 


16.94 


2.74 


1.1120 


17.2 


17.35 


1.03 


1.1108 


17.3 


16.48 


1.07 


1.111 


17.4 


15.97 


1.05 


1.115 


17.5 


16.31 


1.04 


1.111 


17.6 


15.85 


1.09 


1.110 


17.7 


16.04 


1.35 


1.121 


17.8 


17.57 


1.20 


1.116 


17.9 


16.87 


1.00 


1.11 


17.10 


17.31 


1.04 


1.1086 


17.11 


17.01 


1.14 


1.098 


17.12 


17.43 


1.12 


1.105 


17.13 


17.51 


0.98 


1.109 


17.14 


17.32 


0.91 


1.112 


17.15 


16.59 


1.07 


1.104 


17.16 


16.92 


1.19 


1.119 


17.17 


17.42 


1.09 


1.115 


17.18 


17.33 


1.11 


1.115 
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Table 3 — Continued 



ID^ [4.5]^ 


phot-z 


spec-z 


17.19 17.29 


1.06 


1.118 


17.20 16.90 


1.18 


1.107 


17.21 16.38 


1.26 


1.115 


17.22= - 


- 


1.114 


17.23= - 


- 


1.110 


ISCS_J1426.1+3403 (^sp) 


= 1.135 


34.1 16.35 


1.03 


1.1439 


34.2 15.75 


1.04 


1.1301 


34.3 16.01 


1.11 


1.1271 


34.4 15.69 


1.00 


1.1328 


34.5 16.12 


1.18 


1.134 


34.6 16.90 


0.91 


1.13 


34.7 17.24 


0.89 


1.144 


ISCS_J1426.5+3339 (^sp) 


= 1.161 


14.1 16.87 


1.10 


1.16 


14.2 15.96 


1.05 


1.157 


14.3 16.44 


1.19 


1.1631 


14.4 16.64 


1.17 


1.1634 


14.5 16.88 


1.04 


1.1637 


ISCS_J1434.5+3427 (^sp) 


= 1.243 '^ 


342.1 16.76 


1.21 


1.240 


342.2 17.71 


1.21 


1.251 


342.3= - 


- 


1.256 


ISCS_J1429.3+3437 (^sp) 


= 1.258 



30.1 


16.65 


1.18 


1.245 


30.2 


16.39 


1.13 


1.26 


30.3 


17.24 


1.15 


1.2576 


30.4 


17.37 


1.07 


1.263 


30.5 


17.41 


1.29 


1.2611 


30.6 


16.86 


1.08 


1.2583 


30.7 


16.97 


1.31 


1.2582 
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Table 3 — Continued 


ID^ [4.5]*^ 


phot-2; 


spec-z 


30.8 15.19 


0.97 


1.2632 


30.9^ - 


- 


1.2546 


ISCS_J1432.6+3436 (^sp) = 1-347 


29.1 17.63 


1.26 


1.3559 


29.2 17.61 


1.25 


1.35 


29.3 16.39 


1.32 


1.35 


29.4 16.06 


1.21 


1.3320 


29.5 17.55 


1.40 


1.347 


29.6 16.77 


1.29 


1.347 


29.7 16.86 


1.30 


1.34 


29.8 16.14 


1.26 


1.353 


ISCS_J1434.7+3519 (^sp) = 1-373 


25.1 17.52 


3.33 


1.37 


25.2 16.98 


1.53 


1.372 


25.3 16.28 


1.34 


1.37 


25.4 15.68 


1.33 


1.374 


25.5^ - 


- 


1.380 


ISCS_J1438.1+3414 (^sp) = 1-413 ^ 


22.1 15.86 


1.25 


1.411 


22.2 15.79 


1.38 


1.418 


22.3 16.91 


1.38 


1.412 


22.4 16.76 


1.35 


1.414 


22.5^ - 


- 


1.412 



^Coordinates will appear in pub- 
lished version. 

^Vega magnitude at 4.5/im; mag 
= 179.5 Jy. 

^Serendipitous spectroscopic 
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source not included in 5a 4.5/im 
catalog. 

'^Additional spectroscopic mem- 
bers for ISCS J1434.5+3427 pub- 
lished in iBrodwin et al.l (120061 ) 



^Additional spectroscopic mem- 
bers for ISCS J1438.1 - f-3414 pub- 
lished in [Stanford et al.l (120051 ) 



